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ABSTRACT 

We combine the latest observations of disk galaxy photometry and rotation curves at moder- 
ate redshift from the FORS Deep Field (FDF) with simple models of chemical enrichment. 
Our method describes the buildup of the stellar component through infall of gas and allows 
for gas and metal outflows. In this framework, we keep a minimum number of constraints 
and we search a large volume of parameter space, looking for the models which best repro- 
duce the photometric observations in the observed redshift range (0.5 < z < 1). We find 
the star formation efficiency to correlate well with Umax so that massive disks are more effi- 
cient in the formation of stars and have a smaller spread in stellar ages. This trend presents a 
break at around i'max ^ 140 km s^^. Galaxies on either side of this threshold have signif- 
icantly different age distributions. This break has been already suggested by several authors 
in connection with the contribution from either gravitational instabilities or supernova-driven 
turbulence to star formation. The gas infall timescale and gas outflows also present a cor- 
relation with galaxy mass, so that massive disks have shorter infall timescales and smaller 
outflow fractions. The model presented in this paper suggests massive disks have formation 
histories resembling those of early-type galaxies, with highly efficient and short-lived bursts, 
in contrast with low-mass disks, which have a more extended star formation history. The ages 
correlate well with galaxy mass or luminosity, and the predicted gas-phase metallicities are 
consistent with the observations of local and moderate redshift galaxies. One option to explain 
the observed shallow slope of the Tully-Fisher relation at intermediate redshift could be small 
episodes of star formation in low-mass disks. 
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1 INTRODUCTION 

Quantifying the star formation history of galaxies constitutes one 
of the major unsolved issues towards a complete understanding of 
galaxy formation. Within the current paradigm of structure forma- 
tion, disk galaxies represent the building blocks from which the 
whole zoo of galaxies is formed. A detailed analysis of galaxy for- 
mation histories is hindered by the complexity of star formation. 
By coupling simple assumptions about star formation and chemical 
enrichment to population synthesis models which track the evolu- 
tion of well-defined stellar populations, one can constrain the possi- 
ble scenarios of star formation in a semi-quantitative way. Further- 
more, by exploring galaxies at high redshift one takes advantage 
of the lookback times which allow us to probe back into the past 
histories of these galaxies. 

Disk galaxies have been complicated to analyze, mostly be- 
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cause of their ongoing star formation, dust and complex geometry 
which prevent us from getting a clear picture of their underlying 
stellar populations. Qualitatively, the photo-spectroscopic proper- 
ties of disk galaxies such as our own Milky Way are compatible 
with a weak and continuous star formation rate (SFR) along with a 
stronger early episode of star formation whose amplitude correlates 
well with galaxy type, being stronger in disks of earlier types (e.g. 
Kennicutt 1998a). 

Studies of disks at moderate and high redshift seem to yield 
a similar formation scenario. Lilly et al. (1998) observed a sam- 
ple of 341 galaxies selected from the CFRS and the LDSS redshift 
surveys, out to z 1, finding bluer U ~ V colours and stronger 
[Oil] emission at moderate redshift, which suggests an increase in 
the SFR by a factor of about 3 from local galaxies out to z ~ 0.7. 
However, the size function was not found to change with redshift, 
and the larger changes are associated with the smaller galaxies 
(however, see Giallongo et al. 2000; Vogt 2000). This is indica- 
tive of earlier formation of massive disks, in agreement with the 
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lower gas fractions found in these systems (e.g. Bell & de Jong 
2000). Boissier & Prantzos (2001) presented a detailed model of 
chemical enrichment in disk galaxies calibrated to the Milky Way 
(Boissier & Prantzos 1999). The application of this model to disk 
galaxies at moderate redshift showed that large disks have already 
completed their evolution by z ~ 1, whereas low-mass disks un- 
dergo a later evolution, resulting in a steepening of the slope of the 
TuUy-Fisher relation with redshift. A similar conclusion is reached 
by Ferreras & Silk (2001) when comparing sets of simple chemical 
enrichment models over a wide range of parameter space using the 
z = optical and NIR colour-UMAx relation as a constraint. The 
star formation efficiency was found to correlate with galaxy mass, 
a point which is questioned by Boissier et al. (2001). However, this 
controversy hinges on the degeneracy between star formation ef- 
ficiency and infall timescale. Both groups agree on the need for a 
rapid buildup of the stellar component in massive disks, a property 
that can be achieved either by a high efficiency of star formation or 
by a shorter infall timescale. 

The evolution of the TuUy-Fisher relation with redshift is a 
powerful discriminator of star formation histories. Vogt (2000) pre- 
sented a sample of 100 faint disk galaxies in the range 0.2 < z < 1 
and no change was found in either shape or slope with respect to 
local samples. A zero point offset < 0.2 — 0.3 mag was found 
in the rest frame B band. Milvang-Jensen et al. (2003) explored 
cluster/field differences in the TF relation at moderate redshift on a 
smaller sample of 8 disk galaxies in cluster MS 1054-03 (z — 0.83) 
as well as 19 field spirals in the z — 0.15 — 0.9 range. At fixed 
rotation velocity, they found cluster spirals to be 0.5 — 1 mag 
brighter. However, their sample is too small to draw a firm conclu- 
sion on the evolution of the Tully-Fisher relation. In a larger sample 
of disk galaxies in 3 clusters at moderate redshift (2 ~ 0.5) Ziegler 
et al. (2003) did not find significant differences with respect to field 
disks. 

The work on 113 disk galaxies in the range 0.1 < z < 1.0 
from the FORS Deep Field on which this work is based (Bohm et 
al. 2004) revealed a striking result, namely a shallow slope of the 
rest frame B band Tully-Fisher relation with respect to local disks. 
Slowly rotating disks (umax ^ 100 km s~^) at moderate redshift 
appear up to 2 magnitudes brighter than their low-redshift counter- 
parts. One should question whether the samples at low and mod- 
erate redshift correspond to a simple one-to-one mapping. Indeed, 
the sample of Lilly et al. (1998) presents a large fraction of bright 
and small disks at z ^ 1 which do not have local counterparts. 

In this paper, we combine photometric data from the FDF disk 
sample with a model which convolves galactic chemical evolution 
histories with stellar population synthesis in order to relate the pho- 
tometric observations with a small set of phenomenological param- 
eters which describe the major drivers of star formation in galax- 
ies. We describe the observations in §2 and our model in §3. The 
results are discussed in §4 and the resulting ages and metallicities 
are presented in §5. Finally, in §6 we discuss our results and give 
the conclusions. Throughout this paper we use a ACDM cosmol- 
ogy with Sim = 0.3 and Ho — 70 km s^^ Mpc^^, as suggested 
by the analysis of the angular power spectrum of the Cosmic Mi- 
crowave Background observed by WMAP (Spergel et al. 2003). 
For this cosmology, the age of the Universe at 2 = and 2 = 1 
corresponds to 13.5 and 5.8 Gyr, respectively. 
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Figure 1. Colour-red.shift plot.s of the sample of FDF moderate redshift disk 
galaxies used in this paper. The filled and empty dots represent galaxies 
with log fMAX > 2.2 and ^ 2.2, respectively. The lines correspond to the 
colours of simple stellar populations from the latest models of Bruzual & 
Chariot (2003), for three ages (from top to bottom): 8, 1, and 0.1 Gyr. For 
each age two metallicities are considered: Zq (solid) and Zq/3 (dashed). 
Notice the bluer colours - which are most sensitive to young stars - suggest 
young stellar ages, whereas the red/NIR colours imply older ages. 

2 OBSERVATIONS 

We use a subsample of the data set presented by Ziegler et al. 
(2002) and Bohm et al. (2004) to study the redshift evolution 
of the Tully-Fisher relation. The galaxies are extracted from the 
FORS Deep Field (FDF), a deep {U BgRIJKs) survey over a 
6' X 6' area (Appenzeller et al. 2000). Spectroscopy with a res- 
olution R ~ 1200 was carried out at FORSl + 2 mounted on the 
VLT. Throughout this paper we refer to the intrinsic rotation veloc- 
ities as Umax since these have been obtained from the flat part of 
the rotation curves. From the complete sample of 77 galaxies with 
measured Umax, we selected those at redshifts 2 ^ 0.5 result- 
ing in a sample of 30 moderate redshift late-type galaxies, which 
cover a range of rotation velocity between 85 and 400 km s~^. 
FigureQshows the colour-redshift diagrams of the 2 > 0.5 sample 
of FDF galaxies used in this paper. The photometry is performed 
over a 2 arcsec aperture, which corresponds to a projected physi- 
cal size of 12 — 16 kpc in the 2 = 0.5 — 1 redshift range. The 
magnitudes were corrected for galactic extinction, following Heidt 
et al. (2003) and intrinsic rest-frame extinction following Tully & 
Fouque (1985), with the convention of a non-negligible extinction 
for face-on disks, so that the minimum value, e.g., in rest-frame 
B-band is Ab ~ 0.27 mag. For more details on the sample re- 
duction and Umax determination see Bohm et al. (2004). In the 
appendix we have explored the effect that a different dust correc- 
tion would have on the model predictions. The lines in figure Q 
show the colours for a set of simple stellar populations from the 
latest Bruzual & Chariot (2003) models for two metallicities: so- 
lar (solid lines) and Zq/3 (dashed lines), assuming three different 
ages (from top to bottom): 8, 1, and 0.1 Gyr. These are predictions 
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Figure 2. A sample of the generic star formation histories used in the mod- 
elling throughout this paper. The SFHs shown here give good fits to the 
colours of a massive {Mg = —22; left panel) and a lower-mass disk 
{Mg = —20; right panel). The top panels show the star formation rate 
{ip{t)), which is driven by a Schmidt-type law fuelled by the infall of pri- 
mordial gas (fit)}. The dashed line (E{t)) shows the contribution to the 
gas component from stars as they reach their endpoints of stellar evolution. 
The lower panels give the evolution of the metallicity with age. The param- 
eters which correspond to each SFH are labelled in the top panels (see text 
for details). The mass rates are normalized to a final stellar mass content of 
10" (left) and 10^° Mq (right). 

for simple populations with a single age and metallicity. Notice that 
the bluer colours - which are very sensitive to young stars - suggest 
young ages, whereas redder colours are more compatible with older 
ages. A more realistic composite model will take into account the 
different contributions from various ages and metallicities to each 
passband. 



3 MODELLING THE SFH OF DISK GALAXIES 

We consider a phenomenological scenario of star formation and 
chemical enrichment based on the models presented in Ferreras & 
Silk (2001). The rationale of our approach is to try to minimize the 
number of parameters describing the buildup of the stellar compo- 
nent. This enables us to run a fairly big number of realizations of 
these models, scanning a large volume of parameter space. We re- 
duce the description of the process of star formation in galaxies to 
a few parameters, namely: 

• Star Formation Efficiency: Star formation is assumed to fol- 
low a variant of the Schmidt law (Schmidt 1959): 

VW = ^Ptot(^)", (1) 
V ptot / 

where pg is the gas volume density and u is a parameter describ- 
ing the star formation efficiency, which is an inverse timescale for 
the processing of gas into stars, ptot is the total baryon density. 



which - for a single zone model - is just a fixed number through- 
out the analysis. The only purpose of including it in equation Q 
is to give u the units of an inverse timescale. Based on the best 
fit to observations from a local sample of normal spiral galaxies 
(Kennicutt 1998b) we decided to use n — 1.5. This exponent 
is also what one obtains for a star formation law that varies lin- 
early with gas density and inversely with the local dynamical time: 

Tp QC pg/tdyn = pgy/Gp^ OC pj'^ 

• Infall: The infall of primordial gas is described by a set of 
two parameters. We assume the infall rate - f{t) - to be a delayed 
exponential function with an infall timescale (r/). We characterize 
the starting point at some epoch by a formation redshift zp, so that 
the age of the Universe when the infall begins is t{zF). Defining 
At = t — t{zF), vje can write the infall rate as: 

/(t)^^Ate-^*/-/ (2) 

This function rises quickly to its maximum value at At — Tf and 
then declines on a slower timescale. After a time Atgo = 3.9r/ 
90% of the total baryon content (ptot) has fallen on the galaxy. 
Notice that this choice of infall rate allows for a very wide of sce- 
narios. Besides the standard model of a quickly rising rate with 
a slowly decaying branch, our parameter space also includes sce- 
narios with an increasing infall rate - i.e. when long timescales 
(Tf ~ 4 — 5 Gyr) and late formation epochs (zp ~ 2) are chosen. 
It is only the data that eventually constrain these scenarios. 

• Outflow fraction: The ejection of gas and metals from winds 
triggered by supernova explosions constitutes another important 
factor contributing to the final metallicity of the galaxy. We de- 
fine a parameter (0^/3^1) which represents the fraction of gas 
and metals from stellar ejecta and lost by the galaxy. Even though 
one could estimate a - model-dependent - value of /? given the star 
formation rate and the potential well given by the total mass of the 
galaxy (e.g. Larson 1974, Arimoto & Yoshii 1987), we leave (3 as 
a free parameter, only to be constrained a posteriori by the data. 

The stellar yields are simplified in the model by a truncated power 
law fitting the results from Van den Hoek & Groenewegen (2002) 
for Intermediate Mass Stars (IMS) (i.e. M < 8Mq). The yields 
from massive stars - which undergo core-collapse - are obtained 
from Thielemann, Nomoto & Hashimoto (1996) assuming solar 
metallicity progenitors. A Salpeter (1955) IMF was used in the 
0.1 - 6OM0 mass range. Other IMFs such as Chabrier (2003) 
were tested, to find no siginificant difference. After all, the slope 
of the IMF in the upper mass region - which mostly controls the 
chemical enrichment - is similar, and the colours do not vary sig- 
nificantly. The most important difference between these two IMFs 
is a systematic change in the stellar mass-to-light ratios. We refer 
the reader to Ferreras & Silk (2001) for details of the model applied 
to disk galaxies. 

Figure |2| shows two sample star formation histories obtained 
by this parametrization. It corresponds to the best fits obtained for 
a typical massive (left) and low-mass disk (right). The parameters 
used for each one are given in the top panels. The star formation 
rate (tp), infall rate (/) and stellar ejecta (E) are shown, normal- 
ized to a final stellar mass content at z — of 10 (left) and 
10^° M© (right). The bottom panels give the evolution of the gas 
phase metallicity. Several models were run with formation redshifts 
between 2 and 5 and the results vary smoothly between these two 
extrema, which correspond to a formation time of 3.2 and 1.2 Gyr, 
respectively. More realistic models should have a well defined cor- 
relation between this formation redshift and the mass of the galaxy. 
However, in our spirit of keeping the number of constraints to a 
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minimum, we ran six grids of models for each galaxy with various 
values of 2f and chose the one with the minimum x^- We want to 
emphasize that the trends presented in this paper are fairly robust 
and insensitive to the particular choice of zp. The remaining three 
parameters {u, (5, tj) are explored over a large volume of parame- 
ter space as described in more detail below. 



3.1 Data constraints 

For each choice of parameters (r/ , (3) we solve the chemical en- 
richment equations which give us a star formation history that is 
subsequently convolved with the latest population synthesis mod- 
els of Bruzual & Chariot (2003) in order to generate a spectral en- 
ergy distribution (sed) at the observed redshift of the galaxy under 
scrutiny (0.5 < z < 1). A function is constructed by compar- 
ing the model prediction with the 2 arcsec aperture photometry. We 
define the function as: 



E 



[c4i)-cT{r)f 



(3) 



where Cz(i) is the observed colour within a 2 arcsec aperture - 
which represents a physical projected distance between 12 and 
16 kpc in the redshift range 0.5 < 2 < 1. The fluxes were ob- 
tained from the coadded FDF images after convolving them to 
the same seeing - see Heidt et al. (2003) for details. The indices 
i = {1, 2, 3, 4, 5, 6} correspond to colours U — B, B — g, g — R, 
R — I, I ~ J and J — Ks, respectively. The photometric un- 
certainties in the colours are crz(z) and are typically in the range 
0.05 — 0.2 mag. The model prediction in each filter is given by 
C^(j). For each model realization we compute the rest-frame B- 
band absolute luminosity at the redshift of the galaxy. These lu- 
minosities are normalized with respect to the observed apparent 
brightness in a passband closest to the rest-frame B band. In our 
sample this corresponds to R for 2 < 0.85 and / for 2 0.85 
(Bohm et al. 2004). The total apparent magnitudes correspond to 
the MAG_AUTO algorithm of the Source Extractor package (Bertin 
& Arnouts 1996). 



3.2 Finding the best fit 

Searching for the choice of parameters which give the best fit ac- 
cording to the x^ function defined in the previous section is a rather 
computer-intensive endeavour. We decided to explore a set of 6 
formation redshifts, namely zf = {2,2.5,3,3.5,4,5}. The re- 
maining three parameters, namely {v, Tf, 13) are searched for each 
galaxy in our sample in a 32 x 32 x 32 grid over a wide range of 
values, namely: 

-1.3 ^ log(^/Gyr-i) ^ +1.5, 

< /? s; 1, 
-1.3 < log(T//Gyr) < +0.7. 

The best fit to the data is determined by the maximization of the 
likelihood function £(2f, i^, /9, T/) oc expl—^x^i^Fyi'jPjTf)], 
where x^i^F, /3, Tj) is described above. With the uncertainties 
inherent to any modelling of galactic chemical enrichment and 
given that the constraints used in this paper are only based on a 
few photometric measurements, we caution the reader to take these 
results as semi-quantitative trends. The strongest value of the model 
predictions lie in the relative difference between parameters among 
galaxies. In the following figures we show the best fits for each 
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Figure 3. Predicted model parameters as a function of rotation velocity 
(left) or absolute B-band luminosity (right) for the complete (z > 0.5) 
FDF moderate redshift sample. The filled and hollow circles represent the 
z < 0.75 and 2 0.75 subsamples, respectively. Typical marginalized 
1(T error bars are shown. The parameters are (from top to bottom) star for- 
mation efficiency (1/), infall timescale (Tf), and outflow fraction (13). The 
dashed lines in the top panels represent the median value of the stai' forma- 
tion efficiency, separating the sample at ?^max ~ 140 km s~^. 



galaxy along with the la confidence level obtained by the likeli- 
hood function. 



4 THE STAR FORMATION HISTORY OF DISK 
GALAXIES 

Figure|3]shows the prediction for the best fit values of the parame- 
ters controlling the star formation history as a function of rotation 
velocity (Jeft) or B-band luminosity (rigtit). The sample is sepa- 
rated into a lower redshift (0.5 ^ 2 < 0.75; filled circles) and a 
higher redshift sub-sample (2 ^ 0.75; hollow circles). Some typ- 
ical (Ict) marginalised error bars are shown. Even though the un- 
certainties are rather large - as expected from any analysis based 
on broadband photometry - one can still hint at a significant cor- 
relation between the star formation efficiency (z^; top panels) and 
rotation velocity. The dashed lines give the median of 1/ separating 
the sample at Umax ~ 140 km s^ ^ . A clear increase in the star for- 
mation efficiency is seen for the more massive disks. Furthermore, 
the position of this "break" is reminiscent of the recent results of 
Kauffmann et al. (2003) on a large sample of 10^ galaxies from the 
SDSS. Their study showed a clear bimodality in the distribution of 
age-sensitive observables such as the 4000A break strength or US a 
B aimer absorption at a stellar mass i\f* ~ 3 X 1O^°M0, which 
corresponds to log Umax ~ 2.2 for a reasonable choice of stellar 
AI/L. The infall timescale (middle panels) appears also to be cor- 
related with rotation velocity, so that low-mass disks seem to have 
longer infall timescales. The data suggests that v and Tf present 
similar correlations, a well-known degeneracy (e.g. Ferreras & Silk 
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Figure 4. Average (top) and RMS (middle) values of the best fit stellar age 
distribution predicted for each galaxy after being evolved to z = 0. The 
vertical dashed line gives the position of the break at t)MAX = 140 km s~ ^ . 
Notice that galaxies with larger rotation velocities than this break present a 
larger scatter towards older ages. The solid line in the top-right panel gives 
the results from Bell & de Jong (2000) along with a rough estimate of the 
scatter (dashed lines). The bottom panels show the formation redshift which 
gives the best fit (see text for details). 



2003) which would require more detailed photo-spectroscopic ob- 
servables to disentangle. Since this analysis is based on broadband 
photometry which is very much affected by the infamous age- 
metallicity degeneracy, we decided to test the robustness of our 
models with respect to a possible systematic effect from the dust 
correction applied. The appendix compares two possible dust cor- 
rections and proves that the bimodality obtained in the star forma- 
tion efficiency is fairly robust to these corrections. 

At this point we should mention the result of Boissier et al. 
(2001) who claim the efficiency to be uncorrelated with galaxy 
mass. However, their definition of efficiency is different from ours: 
e = il^/AIg, where Mg is the total gas mass. We can translate be- 
tween these two definitions of efficiency: 



(4) 



V AItot I 

The total baryon mass in the system is: Mtot = Mg + Af* + Mej , 
where the last term refers to the gas mass ejected from stars out of 
the galaxy. This term can be written as roughly: M^j ^ PRMi,, 
where (3 is the outflow fraction defined in §3 and R is the returned 
fraction, i.e. the fraction of gas ejected from stars (see e.g. Tinsley 
1980). If we write the observed gas fraction as fg = Mg/{Mg + 
Mi,), then we can rewrite equation|4|as: 



{1 + I3R)- PR 



(5) 



The returned fraction is in the range 0.2 — 0.3 for most IMFs, so 
that we can roughly neglect terms 0{PR) and get the scaling only 
with respect to the gas fraction. Therefore, our results are compat- 
ible with the findings of Boissier et al. (2001) as long as the gas 



Figure 5. Best average stellar (top) and gas phase (middle, bottom) metal- 
licities of the sample. The subsciipt o refers to metallicities measured after 
evolving the best-fit SFH to zero redshift, whereas ^ refers to a measure- 
ment at the observed redshift. The data are shown along with 1-cr error 
bars. The solid line in the top-right panel gives the results from Bell & 
de Jong (2000) along with a rough estimate of the scatter (dashed fines). 
The stars in the left panels are oxygen gas abundances from a sample of lo- 
cal disk galaxies (Gamett 2002), The stars and triangles in the bottom-right 
panel give the gas-phase abundance measured in the z > 0.5 subsample 
of moderate redshift galaxies from CADIS (Maier et al. 2004) and DGSS 
(Kobulnicky et al. 2003), respectively. The lines in the middle-right panel 
represent the luminosity-metallicity relation from a sample of ~ 53, 000 
galaxies from SDSS (thick lines; Tremonti et al. 2004) and from the local 
(K92+NFGS-^KISS) sample explored in Kobulnicky et al. (thin lines; 2003 
and references therein). A ±0.1 dex scatter is assumed. 



fraction decreases with galaxy mass. Bell & de Jong (2000) find 
a correlation between the gas fraction and K band absolute lumi- 
nosity, namely fg — 0.8 + 0.14(Mk — 20), which - assuming 
a constant e and using equation |4|- approximately gives a varia- 
tion in log v of +0.5 dex between bright and faint disks, roughly 
the same range as the one shown in figure|3|for subsamples in the 
same redshift range. Brinchmann & Ellis (2000) studied the stellar 
masses and the star formation rates in a sample of 321 field galax- 
ies out to z ^ 1 and they found a clear decrease of the specific star 
formation rate - defined as i/i/iV/* - with stellar mass (Af*). Fur- 
thermore, ^jj/Mi, was found to increase with redshift. We can write 
the specific SFR in terms of the efficiency: 



M-, 



f. 



1.5 
9 



1-/9 



f9 



(6) 



The anticorrelation found between tp/Mt, and stellar mass is 
thereby a signature of the decreasing gas fraction with galaxy mass. 
Notice that in figure|3|the subsample with the higher redshift (hol- 
low circles) has a systematically higher star formation efficiency 
(v), which is to be expected from equation|6|as the gas fraction will 
increase with redshift at a fixed rotation velocity. 

The bottom panels of figure |3| give the gas outflow fraction, 
which weakly hints at an anticorrelation with rotation velocity. 
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Bright disks are redder and faint disks are bluer This colour trend 
can be explained either by a range of metallicities — in which the 
outflow fraction /3 would play the main role — or by a mixed range 
of ages and metallicities - implying a stronger correlation of v or 
Tf with the luminosity of the galaxy. One can only estimate the 
weight of these two factors by analyses such as the one presented 
in this paper However, the use of a few photometric data points 
to constrain the star formation history prevents us from extracting 
an unambiguous conclusion and only allows us to give a semi- 
quantitative result. The current spectroscopic data from the FDF 
sample cannot be used for these purposes. Indeed, the spectra from 
our z > 0.6 galaxies only feature one emission line - [OII]/3727A 
- with a very low signal-to-noise ratio in the continuum. A detailed 
spectroscopic analysis on a sample such as this one will doubtlessly 
give us a clearer prediction of which mechanisms play a key role in 
the buildup of the stellar component in disk galaxies. 



5 AGES AND METALLICITIES 

For every choice of star formation history as the one shown in 
figure |2| one can compute the moments of the age and metallic- 
ity distributions. Figure |4] presents the average (top panels) and 
RMS (middle panels) ages of the stellar component along with a 
few characteristic la error bars. The filled and hollow circles cor- 
respond to the lower (0.5 ^5 z < 0.75) and higher redshift galaxies 
(z ^ 0.75), respectively. Since the ages and metallicities will be 
strongly dependent on the redshift of each galaxy, we decided to 
put all of them at the same redshift. Hence, the values shown in 
figure|4]correspond to the stellar populations expected after evolv- 
ing the predicted star formation histories to z = 0. The bottom 
panel shows the formation redshift which give the best fit to the ob- 
served colours. Higher formation redshifts seem to be favoured by 
massive disks. The position of the break in the analysis of the star 
formation efficiency is shown as a vertical dashed line. One sees a 
clear increase both in the star formation efficiency and the average 
age as a function of rotation velocity. The massive disks seem to 
have a slightly smaller RMS scatter in the age, a result obviously 
generated by their higher star formation efficiency. 

Figure |5] presents the results for the average mass-weighted 
stellar metallicity (top) and the gas phase metallicity (middle) of the 
best fit models evolved to zero redshift, and the gas phase metallic- 
ity at the observed redshift (bottom). Notice that the stellar age and 
metallicity estimates presented in this paper are in rough agreement 
with the results from the data presented by Bell & de Jong (2000). 
Our data does not allow us to pursue an accurate analysis of the 
ages and metallicities as a function of central surface brightness, 
but we find a similar correlation as Bell & de Jong with respect 
to absolute luminosity. The solid line in the top-right panels of fig- 
ures|4|and|5|show their fits as a function of absolute luminosity. The 
sample was originally presented as a function of AIk and we used 
their published B ~ K colours to show it with respect to AIb ■ The 
dashed lines give a rough estimate of the large scatter found in their 
sample. Gas metallicities are computed assuming a solar abundance 
ratio in order to translate our metallicities (Z) into O/H measure- 
ments. We assume a solar value of 12+ log (O/Jf) = 8.7 (AUende 
Prieto, Lambert & Asplund 2001). A clear luminosity-metallicity 
relation can be seen, which is consistent with local observations 
(K92-I-NFGS-I-KISS as referenced in Kobulnicky et al. 2003; thin 
lines in the middle panel) or with the recent analysis of the oxy- 
gen abundances of a large sample (~ 53, 000) of galaxies from 
SDSS (Tremonti et al. 2004; thick lines). The dashed lines represent 



a ±0.1 dex scatter corresponding to typical error measurements. 
The triangles in the bottom-right panel are gas metallicity measure- 
ments from the Groth Strip survey in the range 0.26 < z < 0.82 
(Kobulnicky et al. 2003) and the star symbols at Mb > —20 corre- 
spond to moderate redshift galaxies from CADIS (Maier, Meisen- 
heimer & Hippelein 2004). Our results are also compatible with the 
local estimates of metallicity from the sample of Garnett (2002), 
shown as star symbols in the bottom-left panel with respect to ro- 
tation velocity. 



6 DISCUSSION 

A simple phenomenological description of the star formation his- 
tory allows us to explore a large number of possible histories only 
to be constrained by the observations. In this paper we have used 
the moderate redshift sample of disk galaxies from the FORS Deep 
Field. Comparing the aperture photometry in several passbands be- 
tween data and simple models of chemical enrichment across a 
wide spectral range (U ■ ■ ■ Ks) we have found an increase of the 
star formation efficiency - as defined in equation^- with rotation 
velocity. A clear threshold appears around vmax ~ 140 km s^^ 
so that slowly rotating disks present uniformly low star formation 
efficiencies and younger stellar populations with a larger spread of 
stellar ages. More massive disks show a wide range of efficiencies. 
A similar conclusion was reached by Ferreras & Silk (2001) on 
the analysis of the local sample of UMa cluster disks from Ver- 
heijen (2001). The thorough analysis of observational data across 
a wide spectral range - from the radio to the UV - probing the 
stellar and gas component of a sample of 928 nearby disk galax- 
ies led Gavazzi, Pierini & Boselli (1996) to conclude that, to first 
order, mass is the driving parameter of galaxy evolution. They sug- 
gested that a positive correlation between star formation efficiency 
and galaxy mass should be invoked to explain the observations. 
Furthermore, this trend is consistent with the analysis of Dalcanton 
et al. (2004) who also find a threshold in Umax when exploring 
the thickness of dust lanes in disks. They relate this threshold to 
a tradeoff between the two main causes of turbulence in the ISM, 
namely supernova feedback and gravitational instabilities. The for- 
mer is the dominant factor in low-mass disks and produces a higher 
velocity dispersion in the gas, which results in lower star forma- 
tion efficiencies. In contrast, Dalcanton et al. suggest turbulence 
in massive disks to be mainly caused by gravitational instabilities. 
The lower velocity dispersion implied by these instabilities results 
in an increased efficiency of star formation. Along similar lines. 
Silk (2001) suggested a Schmidt-like star formation law: 

^^npJ^), (7) 

\VSN / 

where Q, is the disk angular velocity, ag is the cold gas velocity 
dispersion and vsn is a characteristic velocity defined from super- 
nova feedback as the specific momentum injected by supemovae 
per unit star formation rate. In this framework, the results from 
figure |3| suggest a major change in the ratio {ag /vsn) at around 
''MAX ~ 140 km s~^ 

This bimodal trend has been explored in a large dataset of 
~ 120,000 galaxies drawn from the SDSS (Kauffmann et al. 
2003). Low-mass galaxies were found to have younger stellar pop- 
ulations and low concentrations, with some of these systems hav- 
ing experienced recent star formation. At stellar masses above 
~ 3 X 10^" Mq, galaxies feature older populations and the high 
concentrations typical of spheroidal systems. In the light of our 
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Figure 6. The Tully-Fisher relation at moderate redshift. The filled and hol- 
low circles coirespond to the z < 0.75 and z ^ 0.75 FDF subsamples, 
respectively. A few characteristic error bars in log Umax ^2 shown. The 
stars and triangles represent the (z ^ 0.4) data from Vogt et al. (1997) and 
the field disks from Milvang-Jensen et al. (2003), respectively. The solid 
line is a linear fit to the FDF data. The fit to the local sample of Pierce & 
Tully (1992) is also shown as a dashed line, with a shaded region represent- 
ing the scatter 



models, this spheroidal-like behaviour of large disks reveals the 
high star formation efficiency as expected in bulges or elliptical 
galaxies. Furthermore, a recent paper (Kannappan 2004) has ex- 
plored the gas-to-star mass ratios using u ~ K colour as a proxy. 
The sample of ~ 35, 000 galaxies was extracted from the SDSS 
and 2MASS databases and a similar bimodal trend was found, with 
a threshold around a similar stellar mass, i.e. 2 — 3 x lO^^A/Q.Our 
results confirm this bimodal trend in disk galaxies out to z ~ 1. 

The star formation histories proposed in this paper for disk 
galaxies allow us to determine reliable evolution corrections and 
rest-frame luminosities. Figure |6| shows the Tully-Fisher relation 
with respect to rest-frame Mb- The model absolute luminosities 
are very similar to those already presented by Ziegler et al. (2002) 
and Bohm et al. (2004) for which the k corrections were obtained 
using the population synthesis models of Moller et al. (2001). As a 
comparison, we also show (at z ^ 0.4) the field sample of Milvang- 
Jensen et al. (2003; triangles) as well as the sample from the DEEP 
Groth Strip Survey (Vogt et al. 1997; stars). The solid line is a fit 
to our data, whereas the local Tully-Fisher relation is shown as 
a dashed line and the scatter is represented by the shaded region 
(Pierce & Tully 1992). The FDF sample clearly gives a shallower 
Tully-Fisher relation with respect to the local sample. The small 
size and scatter in the sample of Milvang-Jensen et al. (2003) pre- 
vents us from making any estimate of their mismatch with respect 
to the local TFR. Vogt et al. (1997) give a result which is more 
compatible with the observations at 2 = 0. However, their selec- 
tion favours large disks which biases the sample towards early-type 
systems. Bohm et al. (2004) suggested this slope change could be 
caused by a mass-dependent luminosity evolution. In this paper we 



do not attempt to relate the bright and slowly rotating disks from 
the FDF sample with local low-mass disk galaxies. In fact, a naive 
one-to-one mapping of local and high redshift Tully-Fisher rela- 
tions would result in an anticorrelation between star formation ef- 
ficiency and galaxy mass. We only use the aperture photometry of 
the FDF disk galaxies in order to constrain their star formation his- 
tories. The rotation velocities are then used to assign a "size" pa- 
rameter to each galaxy. The evolution of the slowly rotating disks 
from our analysis to z — results in a population of bright disks 
which are not seen in local surveys. Hence, our results conclude 
that these objects do not have local counterparts, in agreement with 
Lilly et al. (1998) who find a population of small and bright disks at 
2: ^ 0.7 which do not have counterparts in the local galaxy census. 
A simple estimate of the contribution from young stars can be done 
if we consider the combination of two simple stellar populations, 
an old one which represents the underlying population and a young 
one with a fractional mass contribution fy which corresponds to 
the recent episode of star formation. We can then relate the actual 
luminosity of the galaxy (L) to the luminosity of the underlying old 
stellar component (Lo), namely 

L _ MT-' / fr X To 

Lo MoTo' yi-fYJTY' ^' 

where T = M/L is the mass-to-light ratio, and the subscripts O 
and Y denote the old and young component, respectively. Let us 
consider a small fractional contribution, say fy = 0.1. Assuming 
a metallicity log{Z/ZQ ) = —0.5 for both components, and an age 
of 6 Gyr for the old population - corresponding to a formation red- 
shift 2f ~ 5 for a galaxy at z = 0.75 - this 10% contribution in 
mass from young stars translates into AAIb ~ Mb ~ Mb,o ~ 
2.5 log(L_B, o/is) ~ —1.9 mag for a young component with 
ty = 0.1 Gyr and —0.95 mag for ty = 0.5 Myr. A couple of Gyr 
after the episode has subsided (i.e. to = 8 Gyr and ty = 2 Gyr) 
the same 10% contribution implies AAIb ~ —0.36 mag. Hence, 
small episodes of star formation could account for the presence of 
these bright and slowly rotating disks. We should mention that no 
strong evolution is expected in wmax from mass accretion from 
z = 1 to the present time. In our adopted ACDM cosmology ac- 
cretion is mostly suppressed at low redshift by the growth factor 
in linear theory. Indeed, the disk galaxy models of van den Bosch 
(2002) show that the rotation curves do not change more than ~ 5% 
in this redshift range. 

Nevertheless, this does not invalidate the correlation between 
vmax and the star formation efficiency, but rather, strengthens it. 
The progenitors of the local slowly rotating disks without any re- 
cent star formation will be below the detection threshold of the FDF 
sample, which implies even lower efficiencies should be considered 
at low Umax, in agreement with the analysis of Ferreras & Silk 
(2001). Their prediction of a steepening of the slope of the TFR 
with redshift implies disks with «max ^ 100 km at z ~ 1 
should be fainter than Mb ~ —18, a result which is also consis- 
tent with the models of Boissier & Prantzos (2001). It is intriguing 
that there is evidence for very modest episodes of late star forma- 
tion even in nearby early-type galaxies (Trager et al. 2000). Theory 
strongly suggests that such episodes should increase with redshift. 
This is commonly accepted as plausible in galaxy compact groups 
and protoclusters, where the merger frequency is high. However a 
comparably strong statement in the environs of typical disk galaxies 
is puzzling. Very few dwarfs are found near the Milky Way or M31 
relative to the predictions of the commonly accepted CDM model. 
The usual explanation is that companion dwarfs have been stripped 
of gas and indeed of stars, as suggested by the evidence for tidal 
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streams (Ferguson et al. 2002) and the ongoing disruption of the 
Sagittarius dwarf galaxy (e.g. Majewski et al. 2003). Such episodes 
must have been more frequent in the past if theory - which predicts 
an order of magnitude more dwarfs than are observed - is to be rec- 
onciled with the observations. Dwarf disruption will include both 
compression of gas in the dwarf core as well as gas ejection into the 
gas reservoir of the host galaxy. Both effects will result in transient 
episodes of star formation. These enhanced phases of star forma- 
tion, although individually modest and short-lived, could account 
for as much as half of the current stellar mass. Such minibursts are 
also inferred to have occurred in the Milky Way from the past his- 
tory of star formation inferred from analysis of chromospheric age 
indicators in the local stellar population (Rocha-Pinto et al. 2000). 
The inferred stages of episodic brightening at a lookback time of 
5 Gyr or more could help account for the brightening that we infer 
of small disks in the past. We conclude that more data are badly 
needed at moderate redshift in the Mb ~ —18 range if we are to 
develop a full understanding of the Tully-Fisher relation. 



APPENDIX A: THE EFFECT OF DUST CORRECTIONS 
ON THE MODEL PREDICTIONS 

Dust correction introduces colour variations which can have an im- 
portant effect in the model predictions. In this paper we have used 
the dust corrections presented in Bohm et al. (2004) which fol- 
low TuUy & Fouque (1985), so that the correction depends on the 
inclination. However, a systematic effect could arise if this dust 
correction were dependent on Umax as presented in TuUy et al. 
(1998). These authors find that luminous disks could suffer more 
reddening than low-mass disks, thereby introducing a systematic 
trend towards older populations in massive disks. In this appendix 
we present the results for a single grid of star formation histories 
choosing a formation redshift zf = 5 and two extinction correc- 
tions: the one used in this paper - only dependent on inclination 
- and the one from Tully et al. (1998) - where both the inclina- 
tion and rotation velocity determine the reddening. For the latter 
we also assume no extinction in face-on galaxies. Figure lATI shows 
the model predictions for the star formation efficiency (top), infall 
timescale (middle) and average stellar ages (bottom) as a function 
of I'M AX. The left panels show the change of these parameters with 
respect to the uncertainty, and the right panels give the parameters 
themselves for these two extinction laws. One can see that the pa- 
rameters vary mostly within the la uncertainties, preserving the 
strong bimodal trend in star formation efficiency. No significant 
systematic trends were seen either by Bohm et al. (2004) when 
comparing these two dust correction prescriptions. Hence, within 
the model uncertainties already discussed in the paper, we conclude 
that the trend seen with respect to star formation efficiency or age 
is imaffected by different dust correction prescriptions. 
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Figure Al. Effect of a different dust correction in the model predictions 
for the star formation efficiency (top), infall timescale (middle) and average 
stellar age (bottom) as a function of Umax- A single formation redshift 
was chosen (zp = 5) with two different extinction corrections. The left 
panels show the variation in the model predictions with respect to the la 
uncertainty. The right panels give separately the model predictions for these 
two corrections, namely: Tully & Fouque (1985: solid dots) and Tully et al. 
(1998; hollow dots). 
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